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A complete analysis of the decoherence properties of a Josephson junction qubit is presented. The 
qubit is of the flux type and consists of two large loops forming a gradiometer and one small loop, 
and three Josephson junctions. The contributions to relaxation (Ti) and dephasing (T$) arising 
from two different control circuits, one coupled to the small loop and one coupled to a large loop, is 
computed. We use a complete, quantitative description of the inductances and capacitances of the 
circuit. Including two stray capacitances makes the quantum mechanical modeling of the system 
five dimensional. We develop a general Born-Oppenheimer approximation to reduce the effective 
dimensionality in the calculation to one. We explore Ti and T$ along an optimal line in the space 
of applied fluxes; along this "S line" we see significant and rapidly varying contributions to the 
decoherence parameters, primarily from the circuit coupling to the large loop. 
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I. INTRODUCTION 

Recent years have seen many successes in obtaining 
high-coherence quantum behavior in a variety of flux- 
based Josephson-junction qubits. The devices which 
show good behavior as qubits are fairly complex elec- 
trical circuits, and a detailed theoretical analysis of these 
circuits has proven useful in arriving at optimal designs 
with the best decoherence behavior 0, 0. Since the first 
reports of coherent oscillations in Josephson qubits 0, 
the observed coherence times have increased by a fac- 
tor of about 5000; theory has had a substantial role in 
this large increase (for a theoretical review of Josephson 
qubits, see 0), by suggesting strategies for choosing op- 
timal settings of control parameters for the operation of 
the qubit. 

In this paper, we report the results of a detailed the- 
oretical study of the flux qubit recently reported by our 
group|5j. We will make extensive use of a method of 
analysis introduced by Burkard, Koch, and DiVincenzo 
(BKD)®. BKD introduced a universal method for an- 
alyzing any electrical circuit that can be represented by 
lumped elements. BKD proceeds in several steps: first, 
the Kirchhoff equations are formulated in graph theoretic 
language so that they describe the dynamics of a general 
circuit in terms of a set of independent, canonical coordi- 
nates. Then, one set of terms in these equations of motion 
(the "lossless" part) is seen to be generated by a Hamil- 
tonian describing a massive particle in a potential; the 
number of space dimensions in which the particle moves 
is equal to the number of canonical coordinates in the 
Kirchhof equations. The "lossy" parts of the equations 
of motion are treated by introducing a bath of harmonic 
oscillators, in the style of Caldeira and Leggett0- 

Finally the resulting total Hamiltonian, involving a 
system, a bath, and a system-bath coupling, can be an- 
alyzed by standard means to determine the decoherence 
parameters, Ti and T^, of the first two eigenlevels of the 



system (the "qubit"). T\ is the energy loss rate of the 
qubit, while T^, the "pure dephasing time", is related to 
the experimental parameter T2, the decay time of Ram- 
sey fringes, by XT 1 = |T 1 ~ 1 +T7 1 . Long T x and T times 
are both necessary conditions for quantum computing. 

The results of this paper have revealed significant facts 
about the dependence of T\ and T^ on the control pa- 
rameters of our qubit. The qubit has, to a good degree 
of approximation, a bilateral symmetry across its mid- 
line (see Fig. This symmetry manifests itself in the 
quantum behavior: The quantum structure is effectively 
that of a symmetric double well potential whenever the 
difference of bias fluxes in the two large loops <£> is the 
flux quantum $ = h/2e. (The structure is a "gradiome- 
ter", meaning that, to good approximation, its behavior 
is only a function of the difference of the magnetic flux in 
the two large loops.) We will analyze the decoherence pa- 
rameters arising from the two impedances shown, Z\ and 
Zi- Because Z\ is coupled to the qubit via the "small" 
loop, we refer to the decoherence parameters associated 
with it as T\ s and T^; the corresponding parameters for 
Z 2 , coupled via the "large" loop, are Tu and T^. We 
find that the bilateral symmetry completely controls the 
overall structure of the XiS and T^s. All these param- 
eters are symmetric in $ around $0 (T^/ and Tu are 
approximately symmetric for small values of control flux 
<f> < 0.39$Oj the other two are exactly symmetric). Fur- 
thermore, Tin, T^s, and T$i all have divergent behavior 
at the symmetric point; T\ s is exactly divergent, T^ s and 
T^i are very nearly so for a large range of small-loop con- 
trol flux $ c . These facts give a powerful motivation for 
operating the qubit always very near <!> = $0 • As a func- 
tion of <& c , T^i is strongly increasing and Tu is strongly 
decreasing (in the symmetric situation). This makes it 
essential to stay within a particular window of operating 
parameters. 

As we will discuss in detail, the full dependence of the 
four decoherence parameters on $ and <£> c is complex, but 
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FIG. 1: The IBM qubit, drawn as a network graph. This 
is an oriented graph of the gradiometer structure coupled 
to two sources Z\ and Zi via mutual inductances M12 and 
M34. Branches of the graph represent Josephson junctions 
Ji , inductances Li and external impedances Zi . Each Joseph- 
son junction branch (thick line) is modelled by a resistively 
shunted Josephson junction (RSJ) containing (inset (a)) an 
ideal junction with critical current I Ct i, junction capacitance 
d and shunt resistance Ri. C4 and C5 denote stray capac- 
itances present in the circuit. The qubit is operated chang- 
ing the external fluxes <E> C and $ applied through the small 
and large loops, respectively. Inset (b): The tree chosen for 
the graph. Values used for this qubit: Lzi = Zi{u))/iuj; 
Z x = Z 2 = 50 n, I c = 2fry = 0.8/xA, C> = lOfF, C 4 = 
C5 = 50fF and, using the modelling program FastHenry, 
L c = 106.27, Li = 32.18, L 3 = L 5 = 605.03, L 6 = L 7 = 
81.46, L 2 = 32.18, L 4 = 605.03, M 12 = 0.8, M 13 = M 15 = 
0.18, Mis = Mvr = -2.41, M 3c = -M 5c = 0.52, M 34 = 
0.5, M 35 = 3.4, M 36 = M 57 = 0.22, Af 37 = M 56 = 
-0.86, Mic = (exact), M 6c = -M 7c = 27.63, M 67 = 
— 13.93 (all in units of pH). 



can be grossly understood as being controlled by two dis- 
tinct regimes, the "semiclassical" and the "harmonic" . In 
the semiclassical regime, the effective potential is a dou- 
ble well with a high barrier between, so that quantum 
tunneling is very small. As <1> C is increased, the barrier 
drops, then disappears altogether; then the qubit poten- 
tial enters the "harmonic" regime, where the potential is 
approximately just a single, quadratic well. These two 
extreme cases are relatively simple; decoherence in the 
regime of crossover between these two is rather complex. 

We gain some further insight into these results via a 
technical improvement that we have added to the analysis 
of BKD 6]. A description of this improvement, an appli- 
cation of the Born-Oppenheimer approximation, is an- 
other important component of this paper. This improve- 
ment was motivated by the fact that we wanted to study 
the effect of stray capacitances in the qubit circuit of 
Fig. n The quantum mechanics that this model defines 
is that of a particle in a five-dimensional potential (five 
because there are three junction capacitances and two 
stray capacitances, each defining a degree of freedom). 
A direct evaluation of the Schroedinger equation in five 



dimensions is numerically complex. But we find that, in 
a controlled way, we can organize these five dimensions 
into four coordinate directions that are "fast" (in which 
the potential rises very steeply) and one that is "slow" 
(and has the double- well structure at low <& c ). Then, just 
as in molecular physics y|, the fast coordinates can be 
treated adiabatically, having the effect of modifying the 
effective slow potential energy in the one remaining co- 
ordinate. The resulting one-dimensional quantum theory 
is very easy to analyze numerically, and amenable to a 
qualitative discussion. 

This paper is organized as follows: Sec. II introduces 
the network graph formalism that we use to analyze the 
quantum mechanics of Josephson circuits. We stress two 
innovations that considerably streamline the analysis: a 
capacitance rescaling, and a Born-Oppenheimer approx- 
imation. The Appendix gives more background about 
the theory, with subsection A giving a review, with some 
minor corrections, of the relevant parts of BKD[6(, and 
subsection B highlighting some new results in network 
graph theory. Sec. Ill discusses the details of the nec- 
essary computation that are specific to the gradiometer 
qubit. Sec. IV gives a qualitative discussion of the fea- 
tures of the four decoherence parameters, T\ s , T\i, T$ s , 
and Tqi that we compute. Sec. V reviews a semiclassi- 
cal analysis from BKD 6] that is helpful in understanding 
the overall features of the decoherence parameters. The 
four Sees. VI-IX give an extended discussion of each of 
the four decoherence parameters. Sec. X gives some con- 
clusions. 



II. ANALYSIS: CAPACITANCE RESCALING 
AND BORN OPPENHEIMER APPROXIMATION 

Our analysis follows closely that of BKD 6J. A sum- 
mary of the essentials of this theory is given in Appendix 
A. The result of this theory is, first, a system Hamilto- 
nian, which we begin with here (see Eq. J73J)): 



H S (t) = -Q^C ^ 



(i) 



+ Vm o¥ > + ^V T [(N * #*)(*) + (S * I B )(t)}. (2) 

To perform the Born Oppenheimer approximation, it is 
best to first go to a rescaled coordinate system in which 
the mass (i.e., the capacitance matrix C) is isotropic. 
This is mentioned in BKD 6] , but we present this analysis 
more generally here to set our notation. We make the 
following coordinate transformation 



q = C 1 /2 C -V2 Qc; 

f = C- 1 /2 C l/2^ 



(3) 
(4) 
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c is some standard capacitance; it is convenient to insert 
this arbitrary number so that q and f have the same units 
as Qc and ip, respectively. Note that the commutation 
relations are left unchanged by this coordinate change: 

7r-(<PiQc,j - Qc,j<fi) = ihSij — > ^-{hq 3 - qjU) = iM i3 . 



The Hamiltonian for the rescaled Schroedinger equation 
is 

HsW = ^q T q+(^)t/'(M), (6) 



C/'(f,t) = -^L^cos(c^( C - 1 /2 f)j) 

i 

+ if T ( c C- 1 / 2 M C- 1 / 2 )f + ^f?y/2 c -i/2 ( N * * B )(t) + c^CT^S * I fl )(t)]. 
2 'Pq 



(7) 



For computing decoherence parameters, we take over un- J_ 
changed the golden-rule formulas discussed in BKD[(J 
(see Appendix): 



\m-<p\0)-(l\m-<p\l)\ 



2 A") 



2k B T. 



(9) 



i = 4|<0|m. V |l)| 2 J(c 01 )coth^L., 



(8) For the rescaled coordinates, these are 



± = ^(OlcVWC-^flDlVMcoth-^-, 



1 



Ic^VC-^flO) - (l|c 1 / 2 m T C- 1 /2 f i 1) |2 



J(w) 



2k B T. 



(10) 
(11) 



We will discuss the use of the Born-Oppenheimer ap- 
proximation to evaluate these formulas. What must be 
computed are matrix elements of the form 



J df(vf)<a|f)(f|, 



(12) 



where a, (3 = 0,1, and v is the constant vector 

c l/2 m T C -l/2_ 

As discussed in the introduction, we single out one 
(more than one is also possible) "slow" degree of freedom 
fjl , and take all coordinate directions orthogonal to this 
one, fj, to be "fast". So 



f = {f||,f±}- 



(13) 



The fast coordinates are characterized by the fact that 
the potential [/'(f) increases very rapidly in the fj_- 
direction; we assume that it is a good approximation to 
expand in these directions to second order: 

[/'(f) « F(f||) + X>(f||)fx,i+^Mf||) f x,ifx,j, (14) 

where b can be taken to be a real symmetric matrix. 

In this case, the Born-Oppenheimer approximation 
is made as follows 0: fix the slow coordinate fn , solve 



the remaining (harmonic) Schroedinger equation in fast 
coordinates fj_. The ground state eigenvalue of this 
Schroedinger equation is 



1 - T b" 1 a+ ^TrVb. 

(15) 

Note that this effective potential has nontrivial fy de- 
pendence from its last two terms. The first and second 
terms represent the value of the potential (in the fj_ co- 
ordinates), and the final term is the sum of the zero point 
energies \hw in this multidimensional harmonic well. 

The minimum of the potential in the fj_ coordinates, 
as a function of fi . is 



fr in ( f n) = -^ b_1 ( f ii) a ( f n)- 



(16) 



The ground state wavefunction in the fj_ coordinates is 
a gaussian centered at this point, which we will indicate 
as 

(fjO,f||)=g(f ± -fr m (f||)). (17) 
In the Born-Oppenheimer approximation, the full wave- 
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function is taken to be 

(f|a) = (f,||a)<f x |0,f,|) = (fy|a) 5 (f ± -fr»(f||)). (18) 

Where (fy |a) is the a th eigenstate of the one-dimensional, 
slow-coordinate Schroedinger equation 



2n\ 2 h 2 d 2 



<f|,|a>=A Q <f|||a). (19) 



We return to the matrix elements that are to be com- 
puted, Eq. (|12fl . We separate the integrand into a fast 
and a slow part: 

J df(vf)(a|f)<f|/3>= / , df(v x .f J .+V||f||)(a|f)<f|/3) 
<H\df±vx ■ f±(a|fy)(f|| |/3).g 2 (f± - ff in (f||)) + 
df I, dfjvii f|| (a|f|| > (f|| |/3) g 2 (f x - f? in (f||)) 
= / ( /r,|( a |f||)(f|||/3) J rffxvx • fxff 2 (fx - fr"(f||)) + 
dE||V||f||(a|f||)<f|||)8) J dfx<7 2 (fx-fr"(f||)) 

= / df||(v ± .fr n (f||))(a|f||)<f|||/3) + 

df||V||f||(a|f||)<f|||/3). (20) 



In the last line we use the fact that the gaussian is a 
normalized transverse wavefunction. The final two-term 
expression of Eq. (|2U|) will be used below in the eval- 
uation of the T\ and expressions. The Schroedinger 
equation solutions (Eq. 1)19(1 ) and all the necessary inte- 
grations are performed numerically in Mathcmatica. 



III. RESULTS FOR THE GRADIOMETER 
QUBIT 

We have calculated the coherence properties of the gra- 
diometer qubit of Koch et al. jj], assuming coupling to 
two different lossy circuits, one inductively coupled to the 
small loop, and the other inductively coupled to one of 
the large loops (see Fig. 0. Here we do not include the 
additional structure considered in a low- loss termi- 
nated transmission line inductively coupled to the other 
large loop (not shown). This structure strongly modi- 
fies the quantum behavior of the qubit when the energy 
splitting of the ground and first excited state of the qubit 
is large (comparable to 1.5GHz, a typical resonant fre- 
quency for the terminated transmission line); however, 
for smaller energy gaps this structure is expected to be 
unimportant. The two lossy structures included are ex- 
pected to account for most of the dissipative and deco- 
hering processes seen by the qubit. 

It is known that the decohering effect of two such struc- 
tures is non-additive, see Brito and Burkard (Ref. 




FIG. 2: Contour plot of the potential 17' (f) on the S line for 
the external fluxes $ c = 0.36<£>o and $ = $o- The red dashed 
line indicates the "slow" direction fy . Along this direction the 
potential is a symmetric double well, with the two relevant 
minima of the potential indicated by dots. The bars show the 
spatial extension of the wave function, in the vicinity of the 
minima, in the "fast" direction fx with the smallest curvature 
of the potential. 



but they show that this nonadditive effect is typically 
small, and we will consider the irreversible effects of each 
structure separately. 

We have extended the analysis of Q to include the 
effect of stray capacitances on the qubit quantum behav- 
ior. We approximate the distributed stray capacitances 
as two new lumped circuit elements, shown with dotted 
lines in Fig. ^ Including these capacitors, the circuit 
theory leads to a quantum description of the qubit that 
is equivalent to that of a particle in a five-dimensional po- 
tential. Using the Born-Oppcnhcimcr analysis developed 
in this paper, the complexity of the calculation is not too 
greatly increased by these additional capacitances. As 
we will see, these extra capacitances, even though their 
capacitances are larger than the junction capacitances, 
cause only quantitative differences in the behavior of the 
decoherence parameters. 

Fig. |2 shows a two-dimensional slice of the potential 
U', after rescaling the capacitance matrix as indicated in 
Eqs. J3JHJ). The slice is chosen to include the two eigendi- 
rections of the rescaled curvature matrix of the quadratic 
part of the U' potential (C^MoC" 1 / 2 of Eq. 0. In 
one of these directions the curvature is zero; in this di- 
rection only the Josephson energy is nonzero, and the 
potential is periodic (about two periods are shown in the 
figure) . This periodicity reflects the 27r periodicity of the 
superconducting phase of the central island of the circuit 
(the place where J\ , J2 , and J3 meet in Fig. . The dis- 
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FIG. 3: Plots of the potential U' , in the vicinity of the rel- 
evant minima, for each one of the orthogonal components of 
f. These plots were done by keeping fixed four coordinates 
at their minimum points and varying the other one. For the 
"slow" direction fn (solid line) we see a double well struc- 
ture (symmetric on S line). Along the "fast" coordinates 

fj^, here calculated at the LI point, a an almost harmonic 
well is present. The external fluxes used for these plots were 
$ c = 0.36$ and $ = $ . 

placement of the two-dimensional plane shown in Fig. [21 
is chosen so that the inductive energy is minimized — re- 
call that the inductive energy consists of a quadratic and 
a linear part. The two dots in Fig. [3 indicate the mini- 
mum energy points of the potential in this plane, which 
is almost (but not precisely) the position of the absolute 
minima (these have also a small component in the other 
three coordinate directions). We choose the "slow" co- 
ordinate f|| of the Born-Oppenheimer approximation to 
be along the line connecting the two minima in the plane 
shown; the other four directions are treated as the "fast 
coordinates" . 

Fig. |31 gives more detail about the potential in these 
"fast" directions. As expected, the potential rises more 
steeply in all these directions than in the "slow" direction. 
The potentials are all basically harmonic, with some no- 
ticeable anharmonicity, particularly in the softest "fast" 
direction fj 1 . But a calculation of the extent of the 
ground wavefunction in this direction (error bars near 
L2 and R2 in Fig. |5J) shows that it remains well confined 
within the harmonic region. 

We have chosen a "symmetric" setting for the param- 
eters, $ = such that the potential is a symmetric 
double well — the depth of the pair of potential min- 
ima in Fig. |21 is equal. This defines a line in the <&-$ c 
plane that we refer to as the "S line" (S for symmetric). 
As the external control parameters $ and $ c are varied, 
this potential landscape is changed in two different ways: 

1. As $ c is varied, the distance between the two min- 
ima, and thus the height of the barrier separat- 
ing them, varies. Increasing $ c from the value 



FIG. 4: The energy difference between the ground and first 
eigenstate, huoi = (l|Hs|l) — (0|Wg|0), as a function of the 
"control" flux $ c on S line. Inset: Detailed view. The red 
dashed line represents a exponential fit of the data, giving a 
reasonable representation of the data in this vicinity. 

shown, <& c = 0.36$o, the distance between LI 
and Rl (L and R for "left" and "right") drops 
rapidly, as shown by Figs. 161 101 which show how 
these minimum points evolve as a function of $ c 
along the S line. As the minima approach one an- 
other, the height of the barrier separating the LI 
and Rl minima decreases rapidly, as shown in Fig. 
1111 In this regime the quantum-mechanical tun- 
nel splitting between the lowest-lying energy lev- 
els increases dramatically, see Fig. 0] Around 
$ c = 0.39<i>o the barrier vanishes entirely. There 
follows a long interval of <i> c in which there is only 
a single minimum per period of the potential; when 
$ c increases a little beyond $ c = 0.39<i>o, the po- 
tential becones quite harmonic around its mini- 
mum. 

2. As $ is varied around <E>0j the energies of the two 
minima are shifted with respect to one another. For 
larger excursions of away from $o> one minimum 
becomes unstable, and only one minimum per pe- 
riod remains stable. 

Fig. [SI shows a large region of the <&-<!> c plane, simulat- 
ing a sequence of measurements very much as they are 
done in the experiment: For a sequence of values of <E> C , 
$ is scanned from left to right and back again. Each scan 
(nearly horizontal line) plots the value of <& c , plus a signal 
proportional to the classical circulating current in one of 
the large loops of the qubit. The most prominent feature 
of this sequence of curves is the thin vertical regions in 
which the scans are hysteretic. This essentially plots the 
region in which there is a double minimum in the po- 
tential. The shape of this region reflects the behavior of 
the barrier height with control flux, Fig. ^2 Looking at 
flux $ = $o, one sees, as one decreases the control flux 
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FIG. 5: Simulated scans of the critical current over a wide 
region of the <1>-$ C plane. These scans periodically show hys- 
teresis in vertically oriented regions in this plane, indicating 
the presence of the double minimum of the potential in these 
regions. 



$ c from about 0.4<I>o, a rapid widening of the hysteresis 
feature, reflecting a rapid increase of the barrier height. 
The abrupt switch to shrinkage of the hysteresis loop re- 
flects a switching of the lowest barrier from the Ll-Rl 
line to the L1-R2 line. This cuspy feature is readily seen 
in the experiment [To| . and is an excellent landmark for 
calibrating the actual applied fluxes. Fig. [5] is clearly 
periodic with changes in applied flux. Since $ and <£> c 
are Aharonov-Bohm fluxes (i.e., involving no magnetic 
field penetrating the interior of the conductors), chang- 
ing either by an integer multiple of $o should leave the 
quantum behavior of the system invariant. This is ac- 
tually not the periodicity that is seen in Fig. [S] This 
absence of Aharonov-Bohm periodicity, an apparent vi- 
olation of gauge invariance, is a result of the fact that 
the outer perimeter of the qubit is not interrupted by a 
Josephson junction; because the temperature is very low 
compared with the superconducting energy gap, there is 
a very high barrier to the motion of a flux quantum into 
or out of the device. If this barrier is assumed to be in- 
finite (as it effectively is in our model), the states of the 
device fall into noncommunicating sectors. 

Within these sectors, there remains the periodicity 
with respect to varying the external fluxes seen in the 
figure: we can show that if $ is changed by an integer 
multiple of <f>o, (ki — fc 2 — 2fc3)$o (Each ki is any inte- 
ger), the qubit Hamiltonian is invariant if $ c is simulta- 
neously changed by — {k\ + fc 2 ) — L ^+ L (^2 + ^3)- This 
shift of $ and <!> c are associated with the phase changes 
Afi = 2vrfci, Af 2 = 27rfc 2 , Af 3 = 2irk 3 , A£t = ^( - k x + 

fo + z^te+te)) and Af 5 = -^^(fe+fca). The 
inductance factors in these expressions are approximate: 
they are only true in the limit that all mutual inductances 



FIG. 6: The values of the phase associated with the Josephson 
junction Ji at the minima of the potential, as function of the 
"control" flux <I> C , along the S line. Two consecutive pairs of 
minima (Li-Ri) along the periodic direction are shown. Near 
$ c > 0.39, the double -minimum structure collapses rapidly 
to a single minimum per period; in this regime the distance 
between the single minima + 1) is 2n. 



are zero. The pattern of invariance as described by these 
equations is closely matched in experimental data. 

The construction of the quadratic and linear parts of 
the potential in Eq. J7J) require a graph-theoretic analysis 
of the gradiometer circuit, Fig. ^ An appropriate tree 
for the circuit graph is shown in the inset (b) of Fig. ^ 
Using this, the loop matrices defined in the Appendix, 
Eq. 1)450. can be read off by inspection: 
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For the numerical analysis of decoherence parameters, 
we need values for the physical parameters of the circuit. 
For the C matrix, circuit modeling indicates that we can 
take it to be a diagonal matrix with diagonal elements 
{10, 10, 10, 50, 50} (in units of fF). The lOfF capacitances 
are for the Josephson junctions, the 50fF capacitances 
are the "strays". Although the strays are numerically 
the largest capacitances, they do not affect the results 
qualitatively, because of their positions in the circuit. 
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FIG. 7: The values of the phase associated with the Josephson 
junction J2 at the minima of the potential, as function of the 
"control" flux $ c , along the S line. Two consecutive pairs of 
minima (Li-Ri) along the periodic direction are shown. Near 
$c > 0.39, the double -minimum structure collapses rapidly 
to a single minimum per period; in this regime the distance 
between the single minima + 1) is 2n. 



FIG. 8: The values of the phase associated with the Josephson 
junction J3 at the minima of the potential, as function of the 
"control" flux <I> C , along the S line. Two consecutive pairs of 
minima (Li-Ri) along the periodic direction are shown. Near 
$ c > 0.39, the double -minimum structure collapses rapidly 
to a single minimum per period; in this regime the distance 
between the single minima (i)-(i + 1) is 2n. 



The L matrices are denoted: 



L : 



( L x M13 M 15 Afie M 17 \ 
M13 L 3 M 35 M 36 M 37 1 
M 15 M 35 L 3 M m M 57 
Mi 6 M 36 M m U M 67 

\M 17 M37 M57 M 67 U J 



Llk — 
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Lz 2 



(25) 



The numerical values of these parameters are given in the 
caption of Fig. ^ 

The decoherence parameters involve the temperature, 
which we take as T = 5K. This rather high temperature, 
much larger than the bath temperature of a dilution re- 
frigerator, is an accurate reflection of the effective noise 
temperature of the circuits coupled to the qubit. Future 
experiments are planned which will make this effective 
temperature much lower. 

The formal applied flux vector is, <fr x = 
{$ c , $, $ p , 0, 0}. $ c and $ have been introduced 
previously, and $ p is the flux in the third loop, the pick 
up loop. $ p will always to be taken to be zero in the 
analyses here. 



With these matrices we compute the coefficients Mo 
and N(cj) using the formulas in the Appendix (Eqs. 
RH1 (S(ui) does not occur, as no current sources are 
present in the circuit). The applied fluxes are time- 
dependent in the experiments that we are modelling, so in 
principle we need to retain the full frequency dependence 
of N(w). The presence of a frequency dependence in this 
operator is indicative of a retardation phenomenon: the 
Hamiltonian at time t is not a function only of the ap- 
plied fluxes at time t; rather, because of the lossy ele- 
ments in the circuit, H(t) depends on a convolution of 
& x over times preceding t. We find, however, that the 
range in time of the kernel N(i) in this convolution is 
very short: this time range is set by L,2/R.e(Zx(u) = 0)) 
and L4/Re(Z2(w = 0)). For our parameters, this time 
is no more that 10 psec. In experiments the applied 
fluxes are varied on a time scale greater than lOOpsec. 
For this reason, we ignore this retardation effect in all 
our calculations here, and set N(w) = N(w = 0). 



IV. DISCUSSION OF Ti AND 

Figures lT2lll3lll4l andll5lshow the obtained dissipation 
and decoherence rates obtained for the gradiometer qubit 
in the vicinity of the symmetric line, shown as a function 
of changes in the small- and large-loop bias fluxes ($ p is 
taken to be zero throughout). The dependences of these 
quantities is complex, with variations over a large range 
of values (note that all the plots are logarithmic). We 
can explain all the trends seen in these curves. Several 
key facts determine the overall structure of these curves: 
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FIG. 9: The values of the phase associated with the stray 
capacitance Ci at the minima of the potential, as function 
of the "control" flux $ c , along the S line. (The phase of a 
capacitance is proportional to the time integral of the voltage 
across the capacitor.) For this phase all minima pairs (Li-Ri) 
with the same values. This occurs because of the absence of 
a Josephson energy term dependent on this phase. 

• Many of the curves have a break around $ = (1 ± 
<5)<J> , S « 0.01. This is a consequence of a level 
crossing that occurs near this value of 6: for larger 
| S | the lowest two energy eigenvalues of the qubit 
are both in one energy well. Thus, for \6\ > 0.01 
the system is too unsymmetrical for the two qubit 
states to correspond to the left and right wells, and 
consequently the results in this regime are not of 
great interest to us. 

• For small values of the control flux <I> C < 0.39$o 
the barrier is high, and the wave function weight 
is concentrated near the minima of the two wells. 
In this regime, which was referred to as the "semi- 
classical" regime in BKD, the various curves vary 
in predictable ways as the barrier height and well 
asymmetry are changed, as we will detail shortly. 

• For large values of the control flux $ c > 0.39$o 
the barrier vanishes, and the single remaining well 
rapidly becomes almost exactly harmonic. It is 
straightforward to calculate what happens to T\ 
and in this harmonic limit, and we will see that 
the data in this regime can be understood with ref- 
erence to this limit. 

• The lossy circuit coupled to the small loop respects 
the bilateral symmetry of the gradiometer qubit. 
An exact consequence of this is that T± s and T$ s 
are mathematically symmetric around $ = $q- 

• The lossy circuit coupled to the large loop does not 
respect the bilateral symmetry of the qubit. Con- 
sequently, Tn and T$i are not symmetric, but for 
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FIG. 10: The values of the phase associated with the stray 
capacitance C*s at the minima of the potential, as function 
of the "control" flux & c , along the S line. (The phase of a 
capacitance is proportional to the time integral of the voltage 
across the capacitor.) For this phase all minima pairs (Li-Hi) 
with the same values. This occurs because of the absence of 
a Josephson energy term dependent on this phase. 

several separate reasons (different ones in the semi- 
classical and harmonic regimes) these functions, for 
the most part, are very nearly symmetric. Actually, 
if the Born-Oppenheimer corrections to the deco- 
hcrcncc parameters, derived in Sec. [D] were left 
out, Tn would be exactly symmetric. 

• The s curves (Ti s and T^) are very different from 
the I curves (Tu and T^i). This perhaps surprising 
result is explained by the fact that the s functions 
have exactly no contribution from the longitudi- 
nal term in the matrix elements (first term in Eq. 
H20f) ). The longitudinal term usually dominates the 
transverse term (second term in Eq. I|20l) ) when it 
is present, as it is for the I functions. As we will 
see, this makes the character of these curves very 
different from one another. 



V. REVIEW OF SEMICLASSICAL ANALYSIS 

As in BKD, we assume that the potential f '(f) de- 
scribes a double well with "left" and "right" minima at 

fL = {f||L,fr"(f||L)}, (26) 

f« = {fii fl ,fr"(fii iJ )}. (27) 

Then, the semiclassical approximation amounts to as- 
suming that the left and right single-well ground states 
\L) and \R) centered at ?l,r are localized orbitals, hav- 
ing amplitude that vanishes very rapidly away from these 
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FIG. 11: The value of the potential barrier on the S line as 
a function of <£> c . Both the Ll-Rl barrier height (solid black) 
and the L1-R2 barrier height (dashed red) are shown. The 
height of the barrier separating the Ll-Rl decreases rapidly 
as these minima approach one another. For $ c > 0.39, the 
barrier vanishes entirely between LI and Rl. 

minima. Then the two lowest eigenstates can approxi- 
mately be written as the symmetric and antisymmetric 
combinations of \R) and \L), 




where u 01 = VA 2 + e 2 , e = (L\H 8 \L) - (R\H S \R) is 
the asymmetry of the double well, and A = (L\Hs\R) 
is the tunneling amplitude between the two wells. A 
increases almost exponentially with $ c as expected in a 
WKB picture, see Fig. 21 Since \L) and \R) are localized 
orbitals, we approximate the matrix elements (see Eq. 
CD): 

(L|vf|fl)«0, (L|vf|L) ravfz, (i?|vf \R) » vf«. 

(30) 

From Eqs. (|28|l -130 |) the eigenstate matrix elements arc 

(0|vf|l) « iA v . A f, (31) 
2 woi 

(0|v-f|0)-(l|vf|l) « — v-Af, (32) 

wot 

where Af = — fjj. These formulas will be applied 
in different ways to explain the four quantities in Figs. 

Gao 

In this semiclassical approximation with localized 
states, the relaxation and decoherence times both diverge 
if Af can be made orthogonal to v. For a symmetric dou- 
ble well (e = 0), T -> oo for all Af. 



FIG. 12: The relaxation time Ti s associated with the dissi- 
pation source Z\. The T\ a plots are presented as a function 
of changes in the small- and large-loop fluxes, <& c and $ re- 
spectively. As a consequence of a high potential barrier, up 
to $ c « 0.39 the data can be well described by a semiclassical 
model, see Eq. 13711 . For <E> C > 0.39$o, the behavior is nearly 
that in a harmonic potential. The mathematical symmetry 
seen around $ = $o occurs because of the bilateral symmetry 
of the qubit. The break around 5 = ±0.01 (<& = (1 + 5)$o) 
indicates that for larger 5 the lowest two eigenstates are both 
located in one well 

VI. T u 

We see in Fig. ^]that as <I> C increases, T\ s initially is 
almost constant in e and decreasing exponentially with 
<fr c ; this behavior changes fairly abruptly to one which is 
exponentially increasing in $ c , with a sharp maximum 
at e = 0. 

The initial behavior is explained by the semiclassical 
theory. We must specialize the semiclassical theory to a 
fact that is special to the circuit coupling to the small 
loop: as a consequence of the bilateral symmetry of the 
structure, the "naive" longitudinal contribution to the 
matrix elements vanishes, i.e., 

v,| = 0. (33) 

With this, we can specialize the T\ matrix clement Eq. 
(EU thus: 

|<0|v • f|l)| 2 « 1 ^.[vl • (fx(fjf) - fx(ff ))] 2 - (34) 

4 UJ m 

We find that, again as a consequence of symmetry, the 
function f|_(f||) has a special form: at the symmetric- well 
point, it is an even function of fn (assuming the origin 
is centered at the midpoint between the two wells); in 
addition, this symmetry is broken continuously as e is 
made nonzero. This can be summarized by writing the 
start of the Taylor series for v • fj_ (f n ) : 

v - fj_ (f || ) =af 2 + 6ef|| +.... (35) 
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FIG. 13: The relaxation time Tu associated with the dissipa- 
tion source Zi. The Tu plots are presented as a function of 
changes in the small- and large-loop fluxes, 4> c and <3? respec- 
tively. As with Tis, Tu has two distinct regions, the "semi- 
classical" ($ c £ 0.39$ ) and the "harmonic" ($ c > 0.39$ )- 
The approximate symmetry around $ = $o arises from the 
dominant || contributions to the matrix element Eq. 13111 . 
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FIG. 14: The dephasing time associated with the dissipa- 
tion source Z\. The T$ a plots are presented as a function of 
changes in the small- and large-loop fluxes, 4> c and <3? respec- 
tively. The two-peak structure is understood as a manifesta- 
tion of non-equal weight of the states |0) and |1) in the matrix 
elements Eq. 13211 . The mathematical symmetry seen around 
$ = $o occurs because of the bilateral symmetry of the qubit. 



Plugging this in and using f,j* = — for e = gives 



Since f.i ' is slowly varying with $ and <I> C , we have 



|<0|v.f|l)| 2 «i^x6 2 e 2 (ff) 2 , 

4 0J Q1 



Jls{X ^2 £ 2 - A 9 + 9 



1 

A 2 



(36) 



(37) 



This simple functional form fits the curves in Fig. IT51 
very well for $ c = 0.36 - 0.39$ . 

For larger <E> C the trend of T ls is explained by the ob- 
servation that around $ c = 0.39^o, the barrier disap- 
pears and the single minimum rapidly approaches being 
an ideal harmonic potential. If the potential were exactly 
harmonic, with its minimum-curvature direction pointed 
in the || direction, then T\ would diverge. The exponen- 
tial growth of Tis m this regime reflects this approach 
to harmonicity. At all values of $ c it remains true that 
for <f> = $o, Tis is divergent, and the lineshapes around 
$ = $o reflect this. 



VII. T u 

Tu also has two distinct regions, the "semiclassical" 
and the "harmonic" . In both regions, the longitudi- 
nal contributions to the matrix element dominate. This 
means that symmetry-breaking contributions (for <I> = 
(1 ± <5)*I > o) remain very small in all regimes (this is un- 
true for Tfjfi). 

The semiclassical prediction for Tu is 



r^ocKoivfii) 



irfh(tf-f)] 2 - 



(38) 



Tu oc 1 



A 2 ' 



(39) 



This equation predicts a Tu which is exponentially de- 
creasing overall, with a deep minimum at $ = $o, as seen 
in the figure. When the potential becomes harmonic, 
then Tu should approach a constant almost independent 
of <E>, since the harmonic oscillator wavefunctions are only 
shifted by the force proportional to <f> — $o; the matrix 
element is independent of this force. Tu is seen to slowly 
vary with $ c : the variation that is seen presumably re- 
flects the small increase in the harmonic frequency as <fr c 
increases. 



VIII. T^ 

The semiclassical approximation follows the same de- 
velopment as for Tis, with the result (see (|36|l 'l 



|(0|vf|0) 



(llv-fll) 



T<h s oc 



2J2(rR\2 



J 0l 



X &V(f,f 



A 2 



(40) 



(41) 



This last equation predicts a strongly diverging T^ s with 
not very much A (i.e., <!> c ) dependence, as is seen initially 
in Fig. HI 

There is a fairly rapid departure from the semiclassi- 
cal prediction for Ta, s in that the divergence at $ = <E>o 
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FIG. 15: The dephasing time associated with the dissipa- 
tion source Z^. The T$i plots are presented as a function 
of changes in the small- and large-loop fluxes, <J> C and 3> re- 
spectively. The significant breaking of symmetry around the 
S line occurs due the transverse contributions to the matrix 
elements Eq. 1321 . and is associated with the unsymmetrical 
way that this source couples to the qubit. 

splits (symmetrically, as discussed above) into two which 
rapidly move away from the center. This is explained by 
the fact that, once the wavefunctions become somewhat 
delocalized, the difference between the and 1 matrix 
elements in 1)400 becomes nonzero, pushing the diver- 
gence away from the symmetric point. This difference be- 
comes nonzero because the state (the symmetric state 
at $ = $o) has more amplitude between the two minima 
than the 1 (antisymmetric) state. Thus, when weighted 
by the v • fj_ (f m ) function (recall Eq. (JJSJ), the and 1 
matrix elements are (initially) slightly different. 

In the harmonic limit, T^ s should diverge, just as T\ a 
does. Comparing Figs. 1121141 though, shows that the 
details of this divergence are rather different. It is evident 
that as this limit is approached, T^i is dominated by the 
remaining differences of the and 1 matrix elements just 
discussed, which become nearly $ independent as the 
and 1 wavefunctions become more harmonic. 



IX. T 0i 

T^i is shown in Fig. 1151 In the semiclassical regime 
this should be 

ex |(0|vf|0)-(l|v.f|l)| 2 « 4-[ V| |(ff-ff)] 2 , (42) 

w oi 

So 

e 2 + A 2 

cx g2 . (43) 

For small A this predicts, as seen, an almost 
$ c — independent behavior, with a very weak divergence 
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FIG. 16: lis for the case without stray capacitances (C4 = 
C5 =0). We see the same qualitative behavior compared with 
that observed in Fig. 1121 

at <f> = <I>o. As A increases, the divergence gets stronger 
and T^i begins to increase overall. 

In the harmonic limit, again, T^i should diverge. But 
as the || contributions disappear, eventually the trans- 
verse contributions to the matrix elements begin to be 
important. These explicitly break the symmetry, as can 
easily be seen as the shifting of the divergence point in 
the last few T^i curves. 

However, this asymmetry is unlikely to be noticeable 
experimentally. Recall that the physical and T\ are 
(approximately |9j) given by summing the s and / rates. 
The strong asymmetry in T^i occurs only when its con- 
tribution to the rate is very small, and the symmetric T^ s 
will dominate. 



X. DISCUSSION AND CONCLUSIONS 

We conclude with a discussion of the effect of the pres- 
ence of stray capacitances, and on the overall implication 
of our results on decoherence parameters for experiments 
on the gradiometer qubit. Figs. Hfjl and IT7I show the re- 
sults for T\s and T^i for the gradiometer qubit with zero 
assumed stray capacitances. Graphically, the results are 
apparently only slightly changed. This is somewhat an 
impression created by the log scale; a closer examina- 
tion of the T\ result shows that the presence of stray 
capacitances actually improves the relaxation time (i.e., 
makes it longer) by about a factor of 10 in the double- 
well region, while leaving it more or less unchanged in 
the harmonic, single-well region. 

Figs. 1181 and IT§1 provide some explanation for this ob- 
servation. We see that the double-well potential profile 
is very similar in the two cases, with the well depths 
being virtually the same. However, the presence of the 
strays pushes apart the rescaled distance between the two 
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FIG. 17: T^; for the case without stray capacitances (C4 = 
C5 = 0). We see the same qualitative behavior as that ob- 
served in Fig. 1141 

minima. This diminishes the tunneling between the two 
wells, and, not surprisingly therefore, lengthens the re- 
laxation time to go from one well to the other. In Fig. 
ED we see that this effect persists right up to the point 
where the two minima merge at around $ c = 0.39$o. 

Finally, Fig. 1201 gives perhaps the most experimen- 
tally relevant summary of our results for the realistic 
gradiometer qubit parameters (with stray capacitances). 
During qubit operation, it is envisioned Q that the 
qubit will be initialized at small control flux, and will 
then be pulsed rapidly up to high control flux; above 
<f> c = 0.38 — 0.39<i>o, we expect the coherence of the 
qubit to be protected by an oscillator stabilization not 
discussed here. 

The preferred initialization point is at a value of $ c 
well below 0.36<&o. We see that here we have the right 
conditions for initialization of the qubit, in that the T\ 
time will be very long - the figure shows it increasing 
exponentially as <i> c is decreased. This is a simple reflec- 
tion of the very large barrier height in this region. T^, 
and therefore T2 are very small in this region, but this 
is not harmful during initialization. When $ c is pulsed 
upwards, quantum dynamics turn on in a region around 
<i> c = 0.36 — 0.37$o, where the 0-1 frequency is increas- 
ing exponentially through the 100MHz range. This is the 
relevant frequency because it is roughly the inverse of the 
anticipated rise time of the pulse in this region, l-10nsec; 
this is referred to as the "portal" region in [jj. 

Thus, the crucial region for the operation of the qubit 
is in the range of <i> c = 0.37 — 0.38<I>o. It is evidently a 
very perilous region for the qubit: T\ is plunging down- 
ward, dropping to around 200nsec, and T2 is increasing 
from its very small values in the initialization region, but 
does not rise far beyond lOOnsec in this region. Since we 
expect p| to pulse through this region in under lOnsec, 
these times are acceptable for qubit operation; but we 
see that the qubit could not function if it were held in 



FIG. 18: Potential U'{t), Eq. [7| along the "slow" coordinate 
f|j for the cases with (black solid line) and without (red dashed 
line) stray capacitances. The same effective capacitance c is 
used for both. It is observed that the height of barrier is 
almost the same for both cases (the difference is less than 
10 _6 %) while the distance between the minima changes ap- 
preciably. Because of these facts one might expect that the 
tunnelling rate should be higher for the case without "stray" 
capacitances, and thus the Ti should be shorter. This is seen 
in our calculations — compare Figs. H3l and Hr)l The external 
fluxes used for these plots are $ c = 0.36$o and $ = $o- 

this region of $ c without any other protective mechanism 
for a long time. Also, we must beware of other effects, 
such as higher effective temperatures or stronger mutual 
inductance couplings, that would make these times even 
worse. We believe that there have been occasions when 
the conditions of the experiment were worse by a factor 
of 10 or more; in this case, the qubit's coherence is not 
likely to survive even a lOnsec traversal of this region of 
parameter space. 

Fig. [3U|is obtained by adding the inverse relation times 
arising from the small-loop and large-loop circuits (we ig- 
nore the small nonadditive effects explored in Q). T\ is, 
though the whole region of interest, dominated by Ty. 
For small $ c , is also dominated by the large- loop cir- 
cuit; however, for $ c > 0.38$o the small-loop dephasing 
becomes dominant. So, we see that the analysis of both 
circuits is experimentally relevant. 

Finally, the inset to Fig. 1201 shows the effects of imper- 
fections in the setting of <£>, which would put the system 
off the S line. We see that, even for departures of lm$ , 
there are noticeable changes in the decoherence parame- 
ters. T\ is actually increased, reflecting the fact that Tu 
has a minimum at the S line. But the system's sensitivity 
to phase fluctuations increases - is smaller off the S 
line. While the changes seen lm<I>o from the S line are 
not dramatic, these results indicate that much larger de- 
partures should be treated with caution, and subject to 
a full analysis. A concluding word: The results reported 
here have been a useful guide to experiment, but they 
have shown their greatest worth when they are part of the 
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FIG. 19: The "slow" coordinate fn minima, for the relevant 
minima (Li-R;), as a function of $ c along the S line. We 
compare the evolution for the cases with (black solid line) 
and without (red dashed line) stray capacitances. 
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FIG. 20: The total relaxation, dephasing and decoherence 
times (Ti, and T2, respectively) along the S line. We can 
see that (Ti) strongly increases (decreases) as a function 
of <E> C . These facts cause there to be a window of desirable 
operating parameters for the qubit. 



iterative design process itself. Thus, even the rather com- 
plete snapshot given here does not do justice to the full 
role that this analysis has had in the process of perfecting 
the gradiometer qubit. The work has already passed on 
to further questions not touched on here, such as the role 
of additional harmonic oscillator circuits in modifying the 
decoherence parameters 0, E3 , and the problems created 
by introducing qubit-qubit coupling. An intimate rela- 
tionship between theory and experiment will continue to 
be crucial in the continuing development of this qubit 
technology. 
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Appendix: Circuit Theory 

We make extensive use here of the systematic analysis 
of flux qubits initiated in BKD. For completeness, we 
review the formalism presented there, both the network 
graph theory and the Caldeira-Leggett analysis. We have 
added some extensions to this theory, which we separate 
out and present in a separate subsection of this appendix. 



A. Review of BKD 

This subsection is a streamlined summary of the results 
presented in jfj. 

An oriented graph Q = {M , B) consists of N nodes 
J\f = {m, . . . , njy} and B branches B — {b\, . . . , 6s}. In 
circuit analysis, a branch bi — (n a u) , Ubu)) represents a 
two-terminal element (resistor, capacitor, inductor, cur- 
rent or voltage source, etc.), connecting its beginning 
node n a u-) to its ending node riu^ . A loop in Q is a con- 
nected subgraph of Q in which all nodes have degree two 
(the degree of a node n € M is the number of branches 
containing n). For each connected subgraph we choose 
a tree %, i.e. a connected subgraph of Qi which contains 
all its nodes and has no loops. The branches that do not 
belong to the tree are called chords. The fundamental 
loops Ti of a subgraph Qi are defined as the set of loops 
in Qi which contain exactly one chord /, € Qi\7l. 

A complete description of the topology of the network 
is provided by the fundamental loop matrix, defined as 

!1, if bj 6 Ti (same direction as fi), 
— 1, if bj G T (direction opposite to fi), 
0, i: /. ; / 

(44) 

where i = 1, . . . F and j = 1, . . . , B. By labeling the 
branches of the graph Q such that the first N—P branches 
belong to the tree T, where P is the number of disjoint 
connected subgraphs of Q, we obtain 

F< L > = (-F T | 1) , (45) 

where F, the loop matrix, is an (N — P) x (B — N + P) 
matrix. 

The state of an electric circuit described by a net- 
work graph can be defined by the branch currents I = 
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(ii , . . . Ig) , where h denotes the electric current flowing 
in branch b,, and the branch voltages V = (Vi, . . . Vg), 
where Vi denotes the voltage drop across the branch bi. 
If we divide the branch currents and voltages into a tree 
and a chord part, 



I = (ItrJch), 
V = (V tr ,V ch ), 



(46) 
(47) 



Then the Kirchhoff laws can be stated very succinctly 
and universally: 



PI ch = -Itr, 

7T^ r - Vr 



F J V tr 



ch 



(48) 
(49) 



Here & x are the external magnetic fluxes threading the 
loops. 

To write the Hamiltonian of the electrical circuit, we 
must further distinguish the different types of electrical 
circuit elements in the graph. We write 



CJ ?CL ?CR * CZ FCB 

FIT 1 T7 1 T7 1 T7 1 

KJ * KL ? KR ? KZ * KB 



(50) 



The sub-matrices Fxy will be called loop sub-matrices. 
The different chord labels are for Josephson junctions 
(J), linear inductors (L), shunt resistors (R) and other 
external impedances (Z), and bias current sources (B). 
Without loss of generality the capacitors (C) can all be 
taken as tree branches. The tree inductors are labeled 
(K). We note here that in our formalism all capacitors 
should be considered to be in parallel with a Josephson 
junction, even if it is one with zero critical current. 

Finally, to fully define the problem, the electrical char- 
acteristics of each branch type should be defined. The 
current- voltage relations for the various types of branches 
are 



I c siny>, 



(51) 



Qc 

v z (w) 



CVc, 
RIr, 

Z(w)IzM, 



T T 



I'LK 



II 
Ik 



(52) 
(53) 
(54) 

(55) 



Here the diagonal matrix I c contains the critical currents 
7 c ,i of the junctions on its diagonal, and siny is the vec- 
tor (sin<^i, sin<^2> • • • ; sinyjjVj)' Eq. I(52|) describes the 
(linear) capacitors (C is the capacitance matrix), and 
the junction shunt resistors are described by Eq. i|53|) 
where R is the (diagonal and real) shunt resistance ma- 
trix. The external impedances are described by the re- 
lation Eq. I|54|l between the Fourier transforms of the 
current and voltage, where Z(w) is the impedance ma- 
trix. The external impedances can also defined in the 
time domain, 



V z (t)= [ Z(t-r)I z (r)dr=(Z*I z )(t), 

J — oo 



where the convolution is defined as 

rt 



(f*g)(*) 



f (t - T)g(r)dr. 



(56) 



(57) 



Causality allows the response function to be nonzero only 
for positive times, Z(i) = for t < 0. In frequency 
space, the replacement w — > u) + ie with e > guarantees 
convergence of the Fourier transform [Tl| 



(58) 



Z{u) = / Z(t)e luJt dt = / Z(t)e luJt dt 



In our formalism it is necessary to distinguish chord from 
tree inductors, so the inductance matrix must be written 
in block form shown in Eq. 1)55(1 . 

With all these definitions, a universal equation of mo- 
tion for the electric circuit reads 



2tt - 2tt - 

C(p = --Ljhimp - M <p - — (N * - — (S * I B )(t) - R. ip -M d *<p. 

<Dn <J>q 







r 



(59) 



This equation as presented is a slight extension of BKD, 
in that <& x and Ib are allowed to be time dependent. If 
they are time independent, then the expressions for the 
coefficients of this equation of motion are: 

The coefficients of this equation are as follows 



M d (uj) = mL z (cj)- 1 m T , 

Mo = N(w - 0)F£ L , 
N(a; = 0)=Fc L LL 1 LL^. 



(60) 
(61) 
(62) 



M d H=mL z ( [ d)- 1 m T , (63) 

S(uj = 0) = F CB - FcLCL^fF^L^FKB, (64) 

m = F C z-F CL (L^) T F£ L L£F K z. (65) 
With the definitions, as given in Q: 

Lz = Lzz — LzlLll^lz, (66) 
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L L M = L LL - L LZ L ZZ M 1 L 



->ZL, 



Lzl — F^zLrFkl, 



Llz — F kl LkFkz, 



L Z z — Lz + F^zLkFkz, 



(67) 
(68) 
(69) 
(70) 



Performing a standard Born-Markov approximation 
for the system dynamics, one obtains predictions for the 
relaxation times of the system: 



1 

T[ 
1 



= 4|(0|m- y|l)| 2 J(w i)coth 
= |(0|m-^0)-(l|m-^l)| 



^01 
2 J M 



(80) 



2fc B T. (81) 



Lll = L + F kl LkFkl 



(71) 



Fcy — F CY + F CL L 1 LlkL k 1 LkFkyi Y— Z,B, 



Fkl — F K l — L k 1 Ll K , 



Ll — L ( 1l + L 1 LlkLk 1 I j kFkl 



Lk — Lk — LlkL x Llk, 



(72) 
(73) 

(74) 



Lk — Lk (Ik — LkFklL 1 LlkL k 1 ) , (75) 



(76) 



We consider the external impedances contributing to 
dccoherence one at a time (for an analysis of non-additive 
effects, see M^w), which determines the quantities 
in Eq. JBSOHJ, then has the form, 



M d (u) = fiK(u)m.m T , 
K{u>) = L^(w), 

|2 



(82) 
(83) 

H = \m\~, (84) 

m = ih/s/Jl = m/|m|, (85) 
where K (t) is a scalar real function, m is the normalized 
vector parallel to m, and *//Z is the length of the vector 
m (/i is the eigenvalue of the rank 1 matrix mrh T ). Also 



L — L — LlkL k 1 Llk- 



(77) 



The last two terms of Eq. I|59|l describe dissipation, 
handled by the Caldeira-Leggett (see 0), to be reviewed 
shortly. The remaining terms are generated by a (time 
dependent) system Hamiltonian: 

Hs(t) = lQlC- 1 Q c +(^]u(cp,t), (78) 



tf(¥>,t) = -5>^ 



COS tpi 



±<p T M <p + §V[(N * * a )(t) + (S * I B )(t)]. (79) 

I 



( fl ) hnKU). 



B. New results 



(86) 



The following expressions were not contained in BKD, 
and are new to this paper. 

The full frequency dependent expressions for S and N 



S( W )=F CB - F cz L- 1 H(FL-L ZL L-iF^+F CL L- 1 LL- 1 (F^ L -L LZ L- z 1 ( w )Fg z ) L K F KB , (87) 



N(w) = FclLl^Ll 1 ^) ~ F cz Lz 1 (^)Lz L Lzi. 



r 



In our previous work we assumed that 1b and G> x were A final result: it is amusing the write out the full ex- 
time independent, so that only the u> — ► limit of these pression for N(o> = 0), from which M is easily con- 
expressions were presented. structed, in terms of the basic input matrices (loop ma- 
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trices F and inductance matrices L): 



It is clear why it is more manageable to write this ex- 
pression in terms of intermediate quantities; one can see 
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1l + L 1 L L k(Lk^Ll K L 



N(c^ = 0) = F CL 

— Lk (Fkl — L k 1 Llk) L : Llk (Lk — LlkL x Llk) ) Lk (Fkl — L K 1 L 
— LlkLk^lk + F^ L ^1k — Lk (Fkl — L k 1 Llk) L : Llk (Lk — Llk^ X Llk) ) Lk (Fkl — L K 1 L 



j lk) 



-1 T T \ 
LK) 



